THE QUANTUM MACMAHON MASTER THEOREM 



STAVROS GAROUFALIDIS, THANG TQ LE, AND DORON ZEILBERGER 

Abstract. We state and prove a quantum-generalization of MacMahon's celebrated Master Theorem, and 
relate it to a quantum-generalization of the boson-fermion correspondence of Physics. 



1. Introduction 

1.1. MacMahon's Master Theorem. In this paper we state and prove a quantum-generalization of 
MacMahon's celebrated Master Theorem, conjectured by the first two authors. Our result was motivated by 
quantum topology. In addition to its potential importance in knot theory and quantum topology (explained 
in brief in the last section), this paper answers George Andrews's long-standing open problem A of finding 
a natural q-analog of MacMahon's Master Theorem. 

Let us recall the original form of MacMahon's Master Theorem and some of its modern interpretations. 

Consider a square matrix A = (ay) of size r with entries in some commutative ring. For 1 < i < r, let 
Xi := ^2j—i aijXj, (where Xi's are commuting variables) and for any vector (mi, ■ ■ . , m r ) of non-negative 
integers let G(mi, . . . , m r ) be the coefficient of x™ 1 x™ r . . . x™ T in Yll=i X™\ MacMahon's Master Theorem 
is the following identity (see |MM| 1: 

oo 

(1) G(mi,...,m r ) = l/det(7- A). 

m± ,m2 1 ...,m r —0 

There are several equivalent reformulations of MacMahon's Master Theorem; see for example |FZj and 
references therein. Let us mention one, of importance to physics. 

Given a matrix A — (ay) of size r with commuting entries which lie in a ring 1Z, and a nonnegative 
integer n, we can consider its symmetric and exterior powers S n (A) and A n (^4), and their traces tr S n (A) 
andtrA ra (A) respectively. Since 



tr S n (A) = G(mi,...,m r ) 

mi^r ...m-r—n 
oo 

det(J-L4) = ^2(-l) n trA n {A)t n , 



n=0 



the following identity 



(2) : : 7-tt - Y tr S n (A)t n 

in 7£[[f]] is equivalent to Q. In Physics J2j is called the boson-fermion correspondence, where bosons (resp. 
fermions) are commuting (resp. skew-commuting) particles corresponding to symmetric (resp. exterior) 
powers. 
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1.2. Quantum algebra, right-quantum matrices and quantum determinants. In r-dimensional 
quantum algebra we have r indeterminate variables Xi ( 1 < i < r), satisfying the commutation relations 

for all 1 < i < j < r. We also consider matrices A = (ay) of r 2 indeterminates ay, 1 < i, j < r, 
that commute with the x^s and such that for any 2 by 2 minor of (ay), consisting of rows i and i', and 
columns j and j' (where 1 < i < %' < r, and 1 < j < < r), writing a := ay, b := ay, c :— a^j, d :— ayy , 
we have the commutation relations: 

(3) ca = gac, (g-commutation of the entries in a column) 

(4) db = qbd, (^-commutation of the entries in a column) 

(5) ad = da + q~ 1 cb — qbc (cross commutation relation). 

We will call such matrices A right-quantum matrices. 

The quantum determinant, (first introduced in |FRT| ) of any (not-necessarily right-quantum) r by r 
matrix B = (&y) may be defined by 

det,(S) := (-9) _inv(7r) ^ 1 i^ 2 2 • • • b nrr , 

TreS r 

where the sum ranges over the set of permutations, S r , of {1, . . . , r}, and for any of its members, 7r, inv(7r) 
denotes the number of pairs 1 < i < j < r for which iti > itj . 

1.3. A g-version of MacMahon's Master Theorem. We are now ready to state our quantum version 
of Macmahon's Master Theorem. 

Theorem 1. (Quantum MacMahon Master Theorem) Fix a right-quantum matrix A of size r. For 
1 < i < r, let Xi :— Y] aijXj, and for any vector (mi, . . . , m r ) of non-negative integers let G(mi, . . . , m r ) 
be the coefficient of ir™ 1 x™ r . . • x™ r in T][ =1 X™* . Let 

Ferm(A) - ^ (-1)1 J ^det q (Aj) 

JC{l,...,r} 

where the summation is over the set of all subsets J of {1, . . . , r}, and Aj is the J by J submatrix of A, and 

oo 

Bos(A) = G(mi, . . . ,m r ). 

mi ,...,m r =0 

Then 

Bos(A) = l/Ferm(A). 

When we specialize to q = 1 , Theorem^recovers Equation @ , which explains why our result is a q- version 
of the MacMahon Master Theorem. For a motivation of Theorem see Section [3] 

The above result is not only interesting from the combinatorial point of view, but it is also a key ingredient 
in a finite noncommutative formula for the colored Jones function of a knot. This will be explained in a 
subsequent publication, |HL| . 

1.4. Computer code. The results of the paper have been verified by computer code, written by the third 
author. Maple programs QuantumMACMAHDN and qMM are available at: jhttp : //www . math . rutgers . edu/ 
^zeilberg/. The former proves rigorously Theorem^fo 1 " an y fixed r. 



1.5. Acknowledgement. The authors wish to thank the anonymous referee who pointed out an error in 
an earlier version of the paper, and Martin Loebl for enlightening conversations. 
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2. Proof 

2.1. Some lemmas on operators. The proof will make crucial use of a calculus of difference operators, 
developed by the third author in |Z1| . This calculus of difference operators predates the more advanced 
calculus of holonomic functions, developed by the third author in |Z2| . 

Difference operators act on discrete functions F, that is functions whose domain is N r . For example, con- 
sider the shift- operators Mi and the multiplication operator Qi which act on a discrete function F(mi, . . . , m r ) 
by 

(MiF)(mi,...,m r ) := F(mi,...,m 1 _i,ra, + l,m 1+ i,...,m r ) 
(Q i F)(m 1) ...,m r ) := q mi F(m x , . . . , m r ). 

It is easily seen that 

MiQ, = qQiMi. 

Abbreviating Qi by q mi , we obtain that: 

(6) Miq mi =q mi+1 Mi .1 /,'/'" </'" M, for i^j. 

Another example is the operator Xi which left multiplies F by ajj. Notice that XjXi = qXiXj for j > i. 
In the proof below, we will denote Xi by In that case, the identity XjXi = qXiXj for j > i holds in the 
quantum algebra, as well as in the algebra of operators. 

Before embarking on the proof, we need the following readily-verified lemmas. 

Lemma 2.1. (commuting Xi with Xj) For 1 < i < j < r, XjXi — qXiXj. 

Lemma 2.2. (commuting Xi with Xj) For each of the a,ij, define the operator Qij acting on expressions P 
involving ay by QijP{aij) := P(qaij). Then, for any 1 < i,j <r, and integer rrii and any expression F 

xT^XjF = [{Q n l Qj 2 l ■ ■ ■ Qj O,,.: ■ • ■ Q^X^F. 

Lemma 2.3. (Column expansion with respect to the last column): Given an r by r matrix iflij) (not 
necessarily quantum) let Ai be the minor of the entry ai r , i.e. the r — 1 by r — 1 matrix obtained by deleting 
the i th row and r th column. Then 

r 

det q (A) = ^(-g) l ~ r (dct q Ai)a ir . 
i=i 

Lemma 2.4. If A is a matrix that satisfies Equation (jSJ and A' denotes a matrix obtained by interchanging 
the i and j columns columns of A, then det q (A') = (— q)~ mv ^det q (A). 

Proof. Suppose first that we interchange two adjacent colums i and j := i + 1. Consider the involution of 
S r that sends a permutation n to tt' — n(ij). Given n £ S r , let (A; n) = (— l)~ lnv ^' a wi i . . . a Wrr denote the 
contribution of tt in det q (A). Then, det g (^4) = ^^{A^tt). Equation (JSJ implies that 

{A-ix) + {A-iT') = {-q){{A'-Tr) + {A'-ir')). 

Summing over all permutations proves the result when j = i + 1. 

Observe that when j = i + 1, the matrix A' is no longer right-quantum since it does not satisfy JSJ. 
However, the proof used only the fact that © holds for the i and i + 1 columns of A. 

Thus, the proof can be iterated wv(ij) times to commute the i and j > i columns of A. The result 
follows. □ 

Lemma 2.5. (Equal columns imply that det q vanishes): Let A be a right-quantum matrix. In the notation 
of Lemma \2.'A for all j ^ r, 

r 

^2(-qy- r (det q A i )a ij = 0. 

i=l 

Proof. If j = r — 1, it is easy to see that g-commutation along the entries in every column of A imply that 
the sum vanishes. 

If j < r — 1 , use Lemma 12.41 to reduce it to the case of j = r — 1. □ 
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Remark 2.6. One can give an alternative proof of Lemmas 12.41 and 12.51 from the trivial 2 by 2 case and, 
by induction using the g-Laplace expansion of a ^-determinant that is completely analogous to the classical 
case. 

2.2. Proof of Theorem ^ The proof is a quantum- adaptation of the "operator-elimination" proof of 
MacMahon's Master Theorem given in |Z1| . Fix a right-quantum matrix A. 
Observe that G{mi, . . . , m r ) is the coefficient of x\ . . . x® in 

r 

H [mi, . . . , m r ; xi, . . . , x r ) := x r ■ ■ ■ x 2 A x x ll^j • 

i=l 

We will think of H as a discrete function, that is as a function of (mi, . . . , m r ) G N r . H takes values in 
the ring of noncommutative Laurrent polynomials in the XiS, with coefficients in the ring generated by the 
entries of A, modulo the ideal given by 0-101. 

Let's see how the shift operators Mi acts on H . By definition, 

MiH(mi, . . . ,m r ;xi, ...,x r ) = 

T -m, ~ m i+i ^-irii-l -m i _ 1 -mi ymi ynii-iymi+1 1 y-nv 

By moving x^ -1 to the front and Xi in front of Xi mi , and using Lemma 12. II and XjX% = qxiXj, we have 
MiH( mi , ...,m r ;x 1 ,...,x r ) = q ^r+rn r ^+---+m i+1 - mi - m2 ---.- mi ^ x -x^ x -m r . . , x -mi^]X 1 mi • • • X™r. 

By moving Xi next to x^ 1 and using Lemma 12.21 this equals to: 

m r +m r _!H hm, + 1 -mi — m 2 TOj_i 1 

[(Qi2 ' ' ' Qir) mi [Qi\Qa ' ' • Qir)" l2 (QiiQi2 1 Qi'i ' " " Qir)™ 3 ' ' ' (Qjl Qj2 ' ' ' Qi.r- 1 )™ r ^] ' 

X-™' ■■■Xi mi X 1 mi ■■■Xr™', 

which is equal to 

m r +m r _iH hmi+i-mi -m 2 m,-i T -1 

[q auXi t 9 o, l2 x 2 t + <7 ai r x r ) a\mi, . . . , m r , x\, . . . , x r ). 

Multiplying out and rearranging, we get that the discrete function ff(mi, . . . , m r ; a^, . . . , x r ) is annihilated 
by the r operators (i = 1, 2, . . . , r) 

i — 1 r 

7?, := ^ _ g -™ J -2m J+1 aww-i-nwoya^ + (Af 4 - oa)^ + ^ -c m<+2roi + 1+ -'' +2m *- 1+ro 'ay^-. 

j=i j=i+i 

Now comes a nice surprise. Let us define bij to be the coefficient of Xj in Vi- For example, for r = 3 we 
have: 

Mi -an -g mi+m2 ai 2 - g mi+2m 2 +m 3ai3 > 

B= | -q- mi - m2 a 2 i M 2 -a 22 -q m2+ma a 23 

_ q - mi -2m 2 -m 3a3l _ q - m2 -m 3a32 M 3 - a 33 

Lemma 2.7. B is a right- quantum matrix. 

Proof. It is easy to see that the entries in each column of B ^-commute. To prove Equation (JSJ , consider the 
following cases for a 2 by 2 submatrix C of B: C contains two, (resp. one, resp. no) diagonal entries of B, 
and prove it case by case, using the fact that the operators Mj and q mj commute with the ay, and satisfy 
the commutation relations □ 

Now we eliminate x\, x 2 , . . . , x r -\ by left-multiplying Vi by the minor of bi r in B = (bij) times (— q) l ~ r , 
for each i = 1,2, ... ,r, and adding them all up. Since B is right-quantum (by Lemma [2.7(1 . Lemma 12.51 
implies that the coefficients of Xi,. . . ,x r -\ all vanish, and det q (B)x r H = 0. After left multiplying by x" 1 
which commutes with the entries in B, we obtain that 

det q (B)H(mi 1 . . . , m r ; xi, . . . , x r ) = 0. 
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Since the entries of B do not contain Xj's, it follows that det q (B) annihilates every coefficient of H, in 
particular its constant term. Taking the constant term yields 

det q (B)G(mi, . . . , m r ) = 0. 

Here comes the next surprise. 
Lemma 2.8. (a) We have: 

det,(B) = (-l) |J| det,(Aj)Mj 
Jc{l,...,r} 

where J = {1, . . . ,r} — J and Mj = l\ jeJ Mj. 
(b) In particular, 

det q (B)\ Ml= ... =Mr=1 = Ferm(A). 
Proof. Let us expand det 9 (-B) as a sum over permutations tt G S r . We have: 

det q (B) = ^(-«)- inv W^ 1 i^ 22 ---6 3rrr 

TT£S r 

r 

= e n^-?)"™^ 6 ^ 

where inv(7r,i) is the number of j > z such that 7ri > Ttj. Now, 6^ = SijMi — qijCLij, where q^ is a monomial 
in the variables q mk , and ]X ?jrji = 1- Moreover, if Tit = i, then for each j with i < j ^ ftj, the exponent of 
in is 2 if 7Tj < z and if Ttj > i. 
Since ]X = 1, we can move the monomials q^ in the left of HiC - 0)~ mv ^7r<»j ancl then cancel them. 
The monomials commute with all entries of the matrix bij, except with the diagonal ones. Commuting q 2mi 
with bu = 5 Vi iMu — q^ida-ii gives: buq m * — q 2mi (5 7Ti iq 2 Mu — 5^,0^,). In other words, it replaces M, by 
q 2 Mi. Thus, we have: 

r 

TTSZSr- t=l 

= e e n(-9) _inv( " , ' i) ^^ 2inv( ^ i) A^n(-9) _inv(7r " ) (- a -^) 

7r£S r JC{l,...,r}i£J i0J 

Now, rearrange the summation. Observe that every permutation n of {1, . . . , r} gives rise to a permutation 7r' 
on the set {1, . . . , r} — Fix(7r), where Fix(7r) is the fixed point set of tt. Moreover, inv(7r', i) = inv(7r, i) — \{j G 
J : j > Using this, part (a) follows. Part (b) follows from part (a) and the definition of Ferm(A). □ 

Hence 

53 (-l) lJl det q (Aj)MjG( mi , . . . , m r ) = 0. 

JC{l,...,r} 

Summing over N r , we get: 

oc 

53 E (~^ lJl <iet q (Aj)MjG(m 1 ,...,m r ) = 0. 

roi,...,m r =0 Jc{l,...,r} 

For a subset J = {fci, . . . , kj} of {1, . . . , r}, we denote by Gj(m/ Cl , . . . , m^. ) the evaluation G(mi, . . . , m r ) 
at = for all i g" J, and we define 

00 

5,/= 53 G(mi, . . . ,m r ). 

m kl ,...,m kj =0 

Using telescoping cancellation, the inclusion-exclusion principle, and Lemma I2.8f b). the above equation 
becomes 

53 (-l)^Fevm(Aj)Sj = 0. 

JC{l,...,r} 
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Using induction (with respect to r), together with 5*0 = 1, we obtain that Ferm(A)S'{ 1 ,..., r } = 1. This 
concludes the proof of the theorem. □ 

3. Some remarks on the boson-fermion correspondence 

Let us give some motivation for Theorem ^ from the point of view of quantum topology. 
For a reference on quantum space and quantum algebra, see |Kal Chapter IV] and [M], 
Recall that a vector (column or row) of r indeterminate entries x%,...,x r lies in r-dimensional quantum 
space A r \° iff its entries satisfy 

for all 1 < i < j < r. 

Recall that a right (resp. left) endomorphism of is a matrix A — (ay) of size r whose entries commute 
with the coordinates Xi of a vector x — {x\, . . . ,x r ) T £ A r ^° and in addition, Ax (resp. x T A) lie in A r ^°. 
Recall also that an endomorphism of A r \° is one that is right and left endomorphism. 

It is easy to see (eg. in |Kal Thm. IV. 3.1]) that A is a right-quantum (i.e., a right-endomorphism) iff for 



every 2 by 2 submatrix ^° of A we have: 

ca = qac, db — qbd, ad = da + q~ 1 cb — qbc. 

Moreover, A is left-quantum iff for every 2 by 2 submatrix of A (as above) we have: 

ba = qab, dc = qcd, ad = da + q~ 1 bc — qcb. 
Finally, A is quantum iff for every 2 by 2 submatrix of A (as above) we have: 
(7) ba = qab, ca — qac, db = qbd, dc = qcd, cb — be, ad — da + q^ 1 cb — qbc. 

The set of quantum matrices A are the points of the r-dimensional quantum algebra M q (r), which is 
defined to be the quotient of the free algebra in noncommuting variables Xij for 1 < i, j, < r, modulo the left 
ideal generated by the commutation relations of Equation J7J) . 

The algebra M q (r) has interesting and important structure. M q (r) is Noetherian, has no zero divisors, and 
in addition, a basis for the underlying vector space is given by the set of sorted monomials {Yli j a^j 3 \ n ij > 0} 
where the product is taken lexicographically; see |Kal Thm IV. 4.1]. An important quotient of M q (r) is the 
quantum group SL q (r) :— M q (r) / (det q — 1), which is a Hopf algebra |Kal Sec. IV. 6] whose representation 
theory gives rise to the quantum group invariants of knots, such as the celebrated Jones polynomial. 

Observing that 

tr S n {A) = ^ G{m 1 ,...,m r ) 

mi+...m r = n 

trA n (A)= = y det q (Aj) 



JC{1,. ..,r}, | J\=n 



Theorem ^ implies that: 

Theorem 2. If A is in M q (r), then 



1 



oo 



Ferm(A) ^ 



tr S n (A) 



n=0 



Since the algebra M q {r) has a vector space basis given by sorted monomials, it should be possible to give 
an alternative proof of the quantum MacMahon Master Theorem using combinatorics on words, as was done 
in |FZ| for several proofs of the MacMahon Master theorem. We hope to return to this alternative point of 
view in the near future. 
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